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The interacting bosons in one-dimensional inversion-symmetric superlattices are investigated from 
the topological aspect. The complete phase diagram is obtained by an atomic-limit analysis and 
quantum Monte Carlo simulations and comprises three kinds of phases: superfluid, persisted charge- 
density-wave and Mott insulators, and emergent insulators in the presence of nearest-neighbor hop¬ 
pings. We find that all emergent insulators are topological, which are characterized by the Berry 
phase TT and a pair of degenerate in-gap boundary states. The mechanism of the topological bosonic 
insulators is qualitatively discussed and the ones with higher fillings can be understood as a |-filling 
topological phase on a background of trivial charge-density-wave or Mott insulators. 
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I. INTRODUCTION 

Topological phases of matter have bulk gap but gap¬ 
less boundary states urn- These studies have been ex¬ 
tended to interacting fermions and bosons. Specially 
much progress has been made for one-dimensional (ID) 
topological bosonic phases. Classification schemes have 
been proposed for such phases It is recognized that 

long-known ID Haldane phase for spin-1 chain is a topo¬ 
logical bosonic phase [8] . While these studies are impor¬ 
tant for our understanding of topological bosonic states, 
it is also highly desirable to study nontrivial bosonic 
phases in experimental accessible systems. 

Motivated by the experimental progress in studying 
topological phases with ultracold atoms [9], exploring 
topologically nontrivial phases of interacting bosons in 
ID superlattices has attracted intensive attentions. It is 
found that by simply replacing free fermions with inter¬ 
acting bosons in ID topological superlattices [inHH], the 
resulting extended Bose-Hubbard models display non¬ 
trivial topological property [IMZl, which are well un¬ 
derstood from the hard-core limit. As the Bose-Hubbard 
model exhibits rich quantum phases due to multiple occu¬ 
pations of bosons on a single site p^ll2Q] , one may expect 
some novel phenomena emerging in topologically nontriv¬ 
ial superlattice systems beyond the hard-core limit. As 
we shall display in this work, surprisingly plentiful phase 
diagram, including the emergence of topological bosonic 
insulators at both fractional and integer fillings and var¬ 
ious phase transitions induced by varying the chemical 
potential and the strength of superlattice potential, are 
found in the simple superlattice system, which is realiz¬ 
able in cold-atom experiments [211ES] • 

We firstly determine the complete phase diagram of 
the system by a combination of atomic-limit analysis 
and quantum Monte Carlo (QMC) simulations. While 


the phase diagram in the atomic limit contains charge- 
density-wave (CDW) and Mott insulating phases with 
various fillings, these insulating phases persist but are 
separated by superfluid phases when hopping terms are 
included. It is interesting that some nontrivial insulating 
phases emerge between two adjacent persisted insulators 
and these emergent insulators are topologically nontriv¬ 
ial, characterized by the Berry phase tt. On the other 
hand, the persisted insulators are topologically trivial 
as they are adiabatically connected to the insulators in 
the atomic limit. The nontrivial topological property of 
the emergent insulators is further confirmed by the pres¬ 
ence of a pair of degenerate in-gap boundary states un¬ 
der open boundary condition (OBC), which leads to the 
splitting of topological plateau in the p — p curve. By 
varying the superlattice strength, one may observe phase 
transitions between the topologically different insulating 
phases. Our results present a complete understanding 
on the topological phase diagram of interacting bosons 
in ID superlattices, which may shed light on the exper¬ 
imental exploration of the predicted exotic topological 
phases and phase transitions. 

II. MODEL OF INTERACTING BOSONS IN ID 
SUPERLATTICES 

We consider the interacting bosons loaded into the op¬ 
tical superlattice with inversion symmetry in the grand 
canonical ensemble, whose basic physics is described by 
the extended Bose-Hubbard model with a superlattice 
potential: 
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where bj (Sj) is the bosonic annihilation (creation) op¬ 
erator, hj = bjbj is the number operator of bosons, the 
hopping amplitude t is set to be the unit of the energy 
{t = 1), Vj^T = Vj represents a superlattice potential 
with the period T, U represents the strength of on-site in¬ 
teractions, and fi is the chemical potential. The inversion 
symmetry further requires Vj = Vr+i-j with j = 1,..., T, 
which can be realized by a simple bichromatic superlat¬ 
tice potential Vj = Acos{27rj/T + J) by tuning the phase 
to 61^2 = 7r(l — 1/T) or 7r(2 — 1/T), with A the strength 
of the potential [22] . Without loss of generality, we shall 
focus our study on the case with T = 3 and (5 = 27r/3 in 
the following calculations. 

In the hard-core limit = oc, Eq. 0 can be mapped 
to the free fermionic model in the superlattice [20| , which 
has been shown that ID topological phases protected by 
inversion symmetry exist at inversion-symmetric point 
61^2 for fillings of 1/T, (T — 1)/T respectively [23] . 
The ID topological phase is characterized by the Berry 
phase, 7 = ^ A{k)dk, with the Berry connection A{k) = 
i{uk\-^\uk) and Uk the occupied Bloch state [24ll26] . Due 
to the protection of the inversion symmetry, the Berry 
phase 7 mod 2 takes two values: tt for a topological phase 
and 0 for a trivial phase. Corresponding to the nontrivial 
Berry phase tt, there appear a pair of degenerate in-gap 
states under OBC, whose distributions are localized near 
the boundaries. 


III. PHASE DIAGRAM OF SOFTCORE CASE 

It is useful to firstly consider the atomic limit t = 0. 
Whether a boson can be added to the jth site with rij 
bosons is determined by the energy difference AE = 
E{nj -hi) — E{nj) = —ji EVjE 2Unj with E{nj) = 
—jirij + VjUj + Urij {rij — 1) the total energy of the bosons 
on the jth site. If AE < 0, the total energy is lowered and 
one more boson can be added to the site. Thus a series of 
lines determined by —ii-\-VjE 2 Unj = 0 separate different 
insulating phases. The phase diagram of the atomic limit 
in the {A/U^ i^/U) plane is shown in Figj^a). The set of 
parallel lines fi/U = A/U + 2 nj determine the occupa¬ 
tion on the middle site of the unit cell and the separated 
regions have gradually increasing occupation numbers. 
The set of parallel lines ji/U = —A/2t/ + 2nj determine 
the occupation on the side sites. There are two kinds of 
insulators in the phase diagram: Mott insulator with a 
uniform density p = ni (or 77 , 2 ) with ni(n 2 ) the number 
of bosons on the side (middle) sites of the unit cell; CDW 
insulator with a density profile reflecting the modulation 
and an average density p = |ni + |n 2 . Different insu¬ 
lators are characterized by the values of structure factor 
S{Q). For the Mott insulators, the structure factor has 
one peak at Q = 0 with S{Q) = p^. In the CDW insu¬ 
lator, a peak develops at Q = ^ with S{Q) = 
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FIG. 1: (Color online) (a) The phase diagram of softcore 
bosons in the atomic limit, (b) The phase diagram of softcore 
bosons in the presence of NN hoppings, in which the shaded 
colored regions represent topological phases and the dotted 
vertical line is a typical cut along which the detail of the 
simulation will be shown in next figure. The phase of the 
superlattice potential is (5 = when the middle (side) site 
of the unit cell has high (low) potential energy. 

Next we turn on nearest-neighbor (NN) hoppings and 
the phase diagram obtained from QMC simulations is 
shown in Figj^b). Compared to the atomic one, the 
phase diagram is considerably modified. Although each 
insulator in the atomic limit persists, the phase bound¬ 
aries are deformed and incommensurate superfluid re¬ 
gions appear between the commensurate insulating re¬ 
gions. Moreover between two adjacent regions separated 
by the line p/U = — A/2t/ -1- 2nj, a insulator with in¬ 
termediate filling p = {pi P 2 ) /2 emerges ( pi,p 2 are 
the fillings of the two regions). Except of the insulator 
with p = 1/3, the emergent insulating regions are con¬ 
nected to the persisted ones with the same fillings, e.g., 
the shaded regions in FigQb) with p = 1 and p = 4/3, 
and are separated from others by superfiuid phases. 

The above phase diagram is obtained by calculating 
the structure factor 5'(Q), 
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FIG. 2: (Color online) The average density, the superfluid 
density and the structure factor as a function of A along the 
cut with fixed A/U = 0.75 in the phase diagram. 


at large A are topological ones. However for p = | insu¬ 
lators, the Berry phase changes its value from tt to 0 at 
the critical point, and the p = | insulators emerging at 
small A are topological ones. The topological property of 
insulators at other fillings can be analyzed similarly. So 
phase transitions between topologically different insula¬ 
tors can be realized by tuning the superlattice strength. 
Moreover these results provide a complete understand¬ 
ing of the topological properties of the insulating phases 
exhibited by interacting bosons in ID superlattices: all 
emergent insulators in the presence of NN hoppings are 
topological, while those persisted CDW or Mott insula¬ 
tors are trivial since they are adiabatically connected to 
the atomic ones. 


and the superfluid density p^ = with W the winding 
number and ^ the inverse temperature m- An insulator 
is characterized by 5'((3) 7 ^ 0 and p^ = 0 , while a super¬ 
fluid phase by S{Q) = 0 and ps 7 ^ 0. To see the details, 
we present the numerical results on the cut with fixed 
A/U = 0.75 in the phase diagram, along which all the 
typical insulating phases are met. As shown in Figthe 
average density p exhibits various plateaus at commensu¬ 
rate fillings, on which the superfluid density ps vanishes. 
So the plateaus correspond to the imcompressive insula¬ 
tors, whose gaps are the widths of the plateaus. Between 
the insulators, the average density increases continuously 
with the chemical potential and the superfluid density is 
finite, implying the system is in superfluid phase. Among 
the insulators, those persisted Mott and CDW insulators 
are distinguished by the values of the structure factors 
(see the Appendix). 


IV. TOPOLOGICAL PROPERTY OF THE 
INSULATORS 

The topological property of the interacting bosonic in¬ 
sulators is characterized by the Berry phase defined with 
the twisted boundary phase 0 [28l [29] , 

7 = (3) 

where 0 takes values from 0 to 27r and t/e is the cor¬ 
responding many-body ground-state wave-function. In 
FigUa), we calculate the Berry phase as a function of the 
strength of the superlattice potential A at different fill¬ 
ings. The p = I {p = |) insulators are in one piece of the 
phase diagram. For the p = ^ insulators, the Berry phase 
is TT and they are topological. However for the p = | insu¬ 
lators, the Berry phase is 0, thus they are trivial. There 
are two pieces for the insulators with p = 1 (p = |), 
which are connected at a critical point. For p = 1 insu¬ 
lators, the Berry phase changes its value from 0 to tt at 
the critical point, implying the p = 1 insulators emerging 
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FIG. 3: (Golor online) (a) The Berry phase as a function of A 
at several typical fillings, (b) The splits of the plateaus under 
OBG at the filling p = |. (c): (upper) The distribution of 
the bosons on two representative points of the split plateaus 
in (b): black square {p = 1.5), red circle {p = 1.625); (lower) 
the difference Spi compared to the bulk sites. The inverse 
temperature of (c) is /3 = 60. In (b) and (c) the parameter 
A/U = 0.75 is used. 


Due to the bulk-boundary correspondence, there 
should appear a pair of degenerate in-gap boundary 
states under OBC in the topological phases. The plateau 
of the topological insulating phase in the p — p curve 
should be altered: below the critical chemical potential 
corresponding to the energy of the in-gap states, none of 
the in-gap states are occupied and the average density 
changes to p' = p — above the critical chemical poten¬ 
tial, both of the in-gap states are occupied and the aver¬ 
age density changes to p" = p It is verified by our 
QMC simulations. As shown in FigJ^b), the plateaus of 
the topological phases split into two pieces with a jump 
at the critical chemical potential and the magnitude of 
the jump is exactly Moreover the jump tends to be 
vertical in the limit of zero temperature, implying the two 
in-gap states are degenerate. To verify the in-gap states 
are boundary ones, we calculate the distribution of the 
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bosons under OBC on two representative points of the 
split plateaus. As shown in FigJ^(c), the distribution on 
the bulk sites are nearly unchanged, thus the filling of the 
in-gap states happens near the boundaries. When none 
of the in-gap states are filled (the lower plateau), there 
is fractional boson less at each boundary compared to 
the bulk sites. However after both are filled (the higher 
plateau), there is fractional boson more. The results 
for the topological phases with other fillings are similar. 
As a contrast, we do not observe the splitting of plateaus 
and boundary states for trivial insulators. Thus our re¬ 
sults give clear evidence that in the topological phases 
a pair of degenerate in-gap states appear and they are 
localized near the boundaries. 

The nature of the topological phases can be qualita¬ 
tively understood from the aspect of their analog to the 
famous topological trimerized model m, whose configu¬ 
ration is shown in Fig. (a). The main character of the 
model is that the bonds connecting different unit cells 
have larger hopping amplitude, while those inside a unit 
cell are equal and smaller. For the topological bosonic 
insulating phase at | filling, the bosons tend to reside 
on the sites in the minima of the superlattice and aver¬ 
agely there is one boson on two adjacent low energy sites 
[see Fig@ (b)]. Since the two sites are identical, the bo¬ 
son can hop freely between them with a larger amplitude 
to gain more kinetic energy. However the hoppings be¬ 
tween other sites are barriered by the potentials and the 
amplitudes are relatively small. Thus the effective hop¬ 
ping amplitudes form a configuration similar to that of 
the topological trimerized model and the system exhibits 
nontrivial topological property. The topological phase at 
higher filling p is similar except that it has a background 
of CDW or Mott insulator with p = pQ — ^e.g., the 
cases with the fillings p = 1, | shown in Fig 13 (c), (d). 
Furthermore the qualitative picture is consistent with the 
values of the structure factor (see the Appendix). So we 
can term the topological phases with higher fillings as 
topological bosonic CDW or Mott insulators depending 
on the background. 



FIG. 4: (Color online) (a) The topological trimerized model 
for comparison’s purpose, in which the solid (dashed) line 
represents strong (weak) bond. Schematic illustration of the 
superlattice potential and the topological ground state at the 
filling: (b) p = \- {c) p = 1; (d) p = |. 

Finally we emphasize that although the system ex¬ 
hibits emergent insulating phases at other Js where there 
are two identical sites in one unit cell, only those emerg¬ 


ing at inversion- symmetric points are topological. In 
Figj^ we calculate the Berry phase as a function of 5 at 
two typical fillings ioi A/U = 2. It shows that the Berry 
phases of the insulators are only quantized to 0 or tt at 
inversion-symmetric points. However when an insulator 
is emergent, its Berry phase has nontrivial value tt. So the 
identified topological bosonic insulators are protected by 
inversion symmetry, which provide concrete models for 
the classifications of ID topological bosonic phases 016]. 



d/ji 


FIG. 5: (Golor online) The Berry phase as a function of the 
phase 6 of the superlattice at the filling: (a) p = 1/3 ; (b) 
p — 1. Here the parameter A/U — 2 and the lattice size 
L = 9 are used. 


V. SUMMARY 

We study the topological phase diagram of interacting 
bosons in ID superlattices with inversion symmetry. The 
complete phase diagram is obtained and the topological 
properties of the identified insulators are determined. It 
is found that the persisted CDW and Mott insulators 
are topological trivial since they are adiabatically con¬ 
nected to the atomic ones, while all emergent ones are 
topologically nontrivial. We present a qualitative mech¬ 
anism for the topological bosonic insulators. This find¬ 
ing is of interests to cold-atom experiments. The studied 
model represents a simple experimentally accessible sys¬ 
tem and the various topological bosonic phases can be 
realized. One may use Bloch oscillations to measure the 
Berry phase [21], and in situ microscopy to detect the 
boundary states [3T4I33] . 
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Appendix A: More results from QMC simulations 

We use QMC simulations to identify different insula¬ 
tors. Among them the trivial CDW and Mott insulators 
are distinguished by the values of the structure factor 
S{Q) at Q = However one should notice the fact 
that the limit values are reached when A is away from 
the critical points, while the values deviate much from the 
atomic ones near the critical points, as shown in Fig|^ 
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FIG. 6: (Color online) The static structure factor S{Q) at 
Q = ^ as a function of A at different fillings. 

For the emergent nontrivial insulators, the values of the 
structure factor support our qualitative physical picture 
of the topological phases with higher fillings. In the pic¬ 
ture, we can approximate: rij ^ -\-Snj^ with the 

boson number of CDW or Mott insulating backgrounds 
and 5nj that of a ^-filling topological phase. Then we 
have S{Q) ^ Sb{Q) + ‘^Sbt{Q) + St{Q), where 

jk 

jk 

SbAQ) = 

jk 

The value of the structure factor in the topological phase 
away from the critical points can be approximated by the 
above formulas. 

For the phase diagram of the softcore case, we take the 
on-site interaction U = St. We also perform QMC simu¬ 
lations with other values of U. As shown in Figj^ some 
insulating phases tend to vanish when U is decreased, 
which are replaced by the superfluid phases. 

We show the QMC simulations under OBC for the 
topological phase with the filling p = | in Sec.IV. The 
splitting of topological plateau reflects the existence of 
the topological boundary states and the results are gen¬ 
eral for all topological phases. In Figj^ we show two 
other cases with the fillings p = 1, and similar results 

are obtained. 



FIG. 7: (Golor online) The average density of bosons as a 
function of the chemical potential for different values of the 
on-site interaction U. Here A/U = 0.75. 
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FIG. 8: (Golor online) The splits of the plateaus under OBG 
at the filling: (a) p = (b) p = 1. The distribution of the 

bosons on two representative points of the split plateaus and 
the difference Spi compared to the bulk sites: (c) the filling 
p = I and the chemical potential p/U = —0.45, —0.35; (d) the 
filling p = 1 and the chemical potential p/U = 0.8625,1.0375. 
In (c) and (d) black square (red circle) represents the case with 
low (high) chemical potential. 
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